Magnetic field control of intersubband polaritons in narrow-gap semiconductors 
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We investigate theoretically the polariton coupling between the light confined in a planar cavity 
and the intersubband transitions of a two-dimensional electron gas confined in semiconductor quan- 
tum wells in the presence of a vertical magnetic field. We show that in heterostructures made of 
non-parabolic semiconductors, the polaritons do not fit a two-level problem, since the cavity photons 
couple to a non-degenerate ensemble of intersubband transitions. As a consequence, the stationary 
polariton eigenstates become very sensitive to the vertical magnetic field, which thus plays the role 
of an external parameter that controls the regime of light-matter interactions. At intermediate field 
strength we predict that the magneto-polaritons have energy dispersions ideally suited to parametric 
amplification. 
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I. INTRODUCTION 

Intersubband (ISB) polaritons are mixed states formed 
by the strong coupling of the light within a microcavity 
and the intersubband transitions of electrons confined 
in a semiconductor quantum well (QW) embedded in 
the cavity. Since the first experimental demonstration 
in 2003i with a GaAs/AlGaAs multiple quantum well 
(MQW) structure, intense research efforts have been de- 
voted to the study of ISB polaritons. With this kind 
of polaritons the light-matter coupling can reach very 
large values^i^ becoming comparable to (or even larger 
than) the bare frequency of the cavity and of the ISB 
excitations. In this ultrastrong coupling regime, inter- 
esting quantum effects appear 4ri Moreover, since the 
coupling strength is proportional to the square root of 
the number of electrons, it can be controlled by electrical 
gatingi^ Beside the observation of the strong coupling 
regime by means of refiectance spectroscopy as in the 
first experiments, and of photovoltaic measurements 
also the electrical injection of cavity polaritons and their 
electroluminescence is being studied with considerable 
effortjii"— Moreover, the coupling of the ISB transition 
with a surface plasmon supported by a metal grating 
has been demonstrated In the effort of reaching larger 
light-matter couplings toward the ultrastrong coupling 
regime, other materials beside GaAs/AlGaAs have been 
considered, like for instance InAs/AlSb MQWs. Also 
the smaller effective mass of InAs with respect to GaAs 
("^inAs/'^GaAs = 0.39) implies a stronger coupling)^ At 
zero magnetic field, the polaritons can be simply and ef- 
fectively described by a two-level problem;^ where the 
first level is the cavity mode with energy iScav, and the 
second level is the ISB transition with energy E21 be- 
tween the first (ground) and the second (excited) sub- 
band; the coupling is quantified by the Rabi frequency 
ilu, where 2hil.ii gives the splitting of the upper and 
lower polariton branches at the resonance -Ecav =-£'21. 



When a magnetic field B is applied along the QW 
growth axis z, neither the energies nor the strength of the 
ISB-cavity coupling are altered; thus, a fortiori, the two- 
level description of the polariton levels remains valid, if 
we still focus on transitions between the Landau levels be- 
longing to different subbands. A theoretical study on the 
possibility of obtaining ultrastrong magneto-polaritons 
couplings exploiting transitions between Landau levels 
in the same subband is reported in Ref. Il6l . Actually, the 
aforementioned insensitivity to a vertical magnetic field 
is exact only for parabolic-band materials. It remains 
a very good approximation for GaAs-based heterostruc- 
tures, since GaAs shows very little non-parabolicity. On 
the contrary, as we show below, in narrow-gap semicon- 
ductors like InAs or InSb, the band non-parabolicity ef- 
fects cannot be disregarded in the calculation of the po- 
laritonic states. 

In this work, we demonstrate that in the non-parabolic 
case the ISB polaritons cannot be simply described in 
terms of two levels. Instead, the cavity photons couple 
to a non-degenerate ensemble of ISB transitions, giving 
rise to a complex evolution of the polariton dispersion 
for increasing B. We shall show that three different cou- 
pling regimes exist as a function of the intensity of the 
magnetic field. To this end, we consider a InAs/AlSb 
MQW heterostructure grown along the z axis. This 
choice is motivated by the experimental observation of 
ISB polaritons in this system;^ as well by the signifi- 
cant band non-parabolicity of InAs. Band parameters 
and band offsets are taken from Ref. [13. We consider a 
cavity with effective thickness so that for the lowest 
mode = tt/L^ holds. The photon energy is given by 

Eca.v ~ (?)-c/^eoo) y^^"^ + fcz, where fc|| is the in-plane k 

vector and £00 is the dielectric constant of the material 
embedded in the cavity. Due to the usual ISB selection 
rules, we consider only light which is TM polarized (i.e. 
with a component of the electric field along the growth 
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FIG. 1. Subbands of a non-parabolic material at B = (not 
to scale). Since the two subbands have different masses and 
thus different curvatures, the transition energy depends on the 
in-plane wavevector. We calculate: £f — Ei ~ 27 meV and 
AEc = AEk=o - AEkp f» f f meV (while = 3f0 meV). 

axis z). 



II. RESULTS AND DISCUSSION 

The non-parabolicity effect can be described in a QW— 
by an effective mass for the in-plane motion m* , depend- 
ing on the subband index n. We set the QW width to 
6.6 nm and calculate that the ISB transition energy at 
zero magnetic field is = 310 meV, and the effec- 
tive masses for the first and second subbands are to* = 
1.68 and = 2.85 mpe, where = 0.026 toq 

is the InAs bulk effective mass. For the simulation we 
choose the cavity mode coupled to uqw = 5 QWs, each 
with an electron density n2D ~ 5 • 10^^ cm^^ in its first 
subband (n = 1; all calculations are performed at zero 
temperature). At i? = 0, the ensemble of ISB transitions 
forms a finite-width continuum. As depicted in Fig. [1] 
this is due to the fact that the two subbands have dif- 
ferent curvatures (since to* ^ TOj). Therefore, the ISB 
transition energy depends on the in-plane wavevector 
A:|| , reaching its minimum value at the Fermi wavevec- 
tor fell — kp. The width of the continuum is given by 

AEc = 7r?i^n2D(l/™i — l/'7^2). 

With the application of a magnetic field B along the 
growth axis, each subband splits into a set of discrete 
Landau levels (LLs) with approximate energies 

^-AB)-K^{j + l)'^. (1) 

where n is the subband index, j = 0, 1, . . . the LL index, 
and i?" is the subband edge energy at _B = 0. In Eq. ([T]) 
we have assumed that the effective mass depends mainly 
on the subband index n and not on the LL index j, which 
is a reasonable assumption as far as the LL separation is 
much smaller than the intersubband transition energy 
(this was checked for all relevant values of the B field). 
Note that, here and in the following, we do not consider 
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FIG. 2. (Color online) Polariton branches at fixed in-plane 
vector fc|| versus the magnetic field B for a InAs/AlSb MQW 
structure embedded in a microcavity (black solid lines). Blue 
dotted-dashed line: bare cavity mode energy Ecav Red 
dashed lines: ISB transition energies AEj{B), plotted only 
in the B range in which the corresponding LLs in the ground 
subband are not empty. For the parameters, see text. 

explicitly the Zeeman spin splitting of the LLs, since the 
ISB transitions are spin-conserving. 

Within the electric dipole approximation, the ISB tran- 
sitions verify Aj = 0. The transition energy AEj [B) = 
i?2,j — Eij for electrons in the j-th LL is then given by 

^E,iB)=El,^{j^'^neB{l,-l,). (2) 

Wc note that in the parabolic case, since m\ ~ ~ 
TO*, we obtain as expected that AEj{B) = E21 docs not 
depend on the magnetic field: therefore all transitions for 
the different LLs are degenerate at the same energy E21 , 
and we can safely apply the same two-level formalism as 
used at B = 0. 

In the non-parabolic case we have instead an ensem- 
ble of transitions at different B-dependent energies. In 
particular, since TO2 > ml, AEj{B) decreases with B for 
all j values, as shown by red dashed lines in Fig. [5] The 
number of active ISB transitions is given by the number 
of filled LLs in the ground subband n = 1, and thus also 
depends on B. In Fig. [2l each AEj transition energy is 
in fact plotted versus B only in the B range for which 
the En=ij level is not empty, i.e. for B < = ^^^^ ■ i 
for j > 1 (while the j = LL is always filled). 

In order to have a significant coupling with more than 
one ISB level, we choose a cavity geometry with a cav- 
ity mode energy Ecav < E^i (blue dotted-dashed line of 
Fig. [2]), so that in absence of coupling the photon en- 
ergy crosses the bare ISB transition energies. In par- 
ticular, in Fig. [5] we choose an effective cavity thickness 
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Lz = 1.5 /xm and we set eoo = 12.32 for the dielectric 
constant of the InAs cavity. The in-plane wavevector is 
fixed to ~ 4.99 /xm~^, corresponding to an angle of 
propagation inside the cavity oi 9 = 67.2° with respect 
to the z axis. 

For the calculation of the polaritons, we note that each 
allowed transition channel j is independent of the others 
and occurs at a different energy Ai?j (for B ^ Q). We 
thus describe the polariton eigenstates (for a given in- 
plane fc|| vector) as a linear combination of the state \a) 
with one photon in the cavity mode and no ISB excita- 
tions, and the set of states \bj) with one ISB excitation 
associated to a given LL j and no photons in the cavity. 
The coupling between the \a) and \bj) states is then given 
by 

% = (3) 

where Uj is the population of the th LL in the ground 
subband (so that ~ ^20)- The frequency Q. is 

calculated in a way similar to the B ~Q casei^ 

_ / e^7i^n2j7iQiy £'21/21 sin^ 

where /21 is the oscillator strength, which has been 
defined and calculated taking into account the non- 
parabolicity as described in Ref. [l9l: /12 = 0.79 TOo/mre 
for our parameters. The calculated polariton branches 
are represented by black solid lines in Fig. [2j Note that, 
for a given value of B, we have included in the calculation 
only the \bj) states originating from non-empty LLs. 

The system parameters have been chosen in order to 
achieve a significant coupling between the cavity mode 
\a) and more than one ISB transition level \bj). This 
can be achieved only if the coupling energy h(l is of the 
order of the typical deviation of the ISB transition en- 
ergies with respect to Ei^i- In fact, if the coupling is 
much larger than the energy separation between the dif- 
ferent ISB transitions, the latter ones behave essentially 
as a single degenerate level for what concerns the cou- 
pling with the cavity photons, and we then recover the 
ordinary two- level regime (not shown). In InAs/AlSb 
heterostructures, the ISB transition energy deviation is 
typically of the order of 10-15 meV at i? sa 10 T (see 
Fig. [2]). For the parameters used in Fig. [2l the coupling 
?ir2 is about 5 meV. We notice moreover that the effects 
of the squared vector potential can be safely neglected 
in our structure, since the correction is of the order— 
of fi^Qp' / E21, which is much smaller than cavity energy 
Ecav (in our case fiQ. ^ Ecav ~ £'21)- 

From Fig. [2] it is apparent that the polariton levels dis- 
play a complex evolution as a function of the magnetic 
field. We distinguish three field regions. For large B 
fields, where only the j = LL is filled, we recover the 
two-level correspondence, valid also for parabolic materi- 
als. As B decreases, however, more LLs start to be filled 




FIG. 3. (Color online) Squared modulus of the "light" part 
of the polariton eigenvectors. Same parameters as in Fig. [2] 

and as a consequence, more ISB transitions couple to 
the light. At B = Bj (j > 1, 2, . . .) a new state appears, 
which however has a zero coupling at this precise value of 
the magnetic field, since the corresponding LL is empty. 
Decreasing B, its population increases (while the pop- 
ulations of lower LLs decrease), so that the coupling is 
spread between the levels. Finally, in the B — > limit, we 
end up with a bare cavity mode coupled to (and placed 
inside) a finite-width continuum of ISB transitions. As 
it is well knowu)^ the resulting eigenstates depend on 
the ratio between the continuum width and the coupling 
strength. Since in our case AEc ~ 2fi(l, two polariton 
states appear near each side of the bare ISB continuum 
limits (see Fig. [5] and discussion below). 

From the above discussion, we see that the magnetic 
field assumes the role of a real external control parame- 
ter, which can be used to tune the regime of light-matter 
interactions. To illustrate more clearly this point, we 
study also the eigenvector components of the polariton 
eigenstates. In Fig.[3]we show the squared modulus of the 
"light" part of the polariton eigenvectors, i.e. the compo- 
nent of the eigenvectors associated to the state \a). The 
magnitude of this component displays the three different 
regimes mentioned above. At large B fields, we clearly 
identify the two polaritons resulting from the strong cou- 
pling between light and the j = ISB transitions. In 
the opposite B regime, i.e. for small B values, we see 
that the light component is mainly concentrated on the 
two extremal polariton branches. All other polaritons 
have a significantly smaller light component, and thus 
this regime resembles a two-level regime. Note however 
that all states have an influence on the overall coupling 
also in the i? — >■ limit, and therefore they cannot be 
disregarded in the calculation of the polariton coupling. 
Finally, there is a third regime for intermediate values 
of the magnetic field: in this case, the light is coupled 
with a discrete set of ISB states, and the resulting po- 
lariton branches have a similar magnitude of the light 
component. 

To discuss more in detail the intermediate regime, we 
focus on a magnetic field i3 = i?2 = 5.17T, which corre- 
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FIG. 4. (Color online) Polariton branches as a function of the 
in-plane wavevector (black solid lines) for B = B2 =5.17 T. 
The blue dotted-dashed line (red dashed lines) represent the 
cavity mode (j — and j = 1 ISB transitions) in absence of 
coupling. The vertical line indicates the fey = fcyij value (see 
text) at which Fig. [2] is calculated. Orange and green arrows: 
sets of entangled states in an optical parametric oscillator 
phenomenon (see text). Same parameters of Fig. [2l 



sponds to a complete filling of the j = and j = 1 LLs. 
Since the two LLs have the same population, Eq. ([3]) im- 
plies that the respective states |6o) and |&i) couple to the 
light with the same strength Ml/y/2. As a consequence of 
this and of the relative energy position of the uncoupled 
levels, the three resulting polaritons states have a simi- 
lar magnitude of the light component, as it can also be 
deduced from Fig. [31 The dispersion of these three polari- 
tons (at B = B2) a& & function of the in-plane wavevector 
k\\ is shown in Fig. lU the uncoupled cavity mode fre- 
quency and intersubband transitions arc also shown with 
blue dotted-dashed and red dashed lines, respectively. 
The vertical line indicates the value of fc|| used in Fig. [2] 

The magnetic field control of the ISB polaritons might 
be observable in an optical experiment. For very strong 
B fields (not discussed here) the cavity mode is ener- 
getically isolated, well above all the ISB levels, so that 
any spectrum probing the light component of the system 
eigenstates should reveal a single intense peak at the bare 
cavity energy. Decreasing to the high fields of Figs. [2] and 
|3] (i? sa 15 T), the optical spectrum is expected to dis- 
play instead two peaks, characteristic of the strong ISB 
(j — 0)-cavity coupling. The spectrum evolves then into 
three peaks of comparable intensities when decreasing 
the field to B fa B2. More peaks are expected to emerge 
when we further decrease B, if the broadening is small 
enough to allow to resolve them; the central peaks should 
decrease in amplitude as B is further decreased, to the 
benefit of the two main lines at i? = 0. 

Let us finally discuss an interesting aspect of the in- 
termediate field region, resulting from the existence of 
multiple polariton lines. In Fig. |4l for fc|| = (vertical 
gray line) the cavity mode lays exactly at mid-distance 
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FIG. 5. (Color online) Polariton branches for B — B2 = 
5.17 T as a function of the angle of the coupled light in the 
substrate. The blue dotted-dashed line (red dashed lines) 
represent the cavity mode {j = and j = 1 ISB transitions) 
in absence of coupling. The orange arrows connect the same 
states as in Fig. [4] Same parameters of Fig. (2] 



from the two j = and j = 1 ISB transitions. More- 
over, since the bare cavity dispersion is to a good ap- 
proximation a linear function of fc|| around it can 
be easily shown that the resulting polariton dispersions 
Em{k\\) (with m = 1,2,3 for the three branches) have 
the following interesting "mirror" property: they fulfill 
Esik^R + k) + Ei{k\\ji — k) ~ 2i?2(fc||_R), with k a small 
deviation from the resonance wavevector. Two of such 
sets of three-polariton states are pictured by green and 
orange arrows in Fig. The three polaritons of each 
set are thus exactly phase-matched in both energy and 
wavevector spaces. Additionally, the upper and lower 
states have always identical group velocities, while all 
three waves velocities coincide for a particular value of 
detuning k = ky (orange arrows in Fig. |4]). This might 
lead to improved non-linear optical response, like in the 
optical parametric oscillation phenomenon;^ which has 
been studied in the literature in monolithic semiconduc- 
tor microcavities exploiting exciton polaritons^"— and 
coupled microcavitiesi^ After pumping on the central 
state, entangled photon pairs (idler and signal) would be 
expected from the upper and lower branches. These lat- 
ter would propagate along well-defined directions with 
respect to the pump beam, allowing an angular discrim- 
ination of the beams at the sample outcome (see below) . 
Moreover, the generation would be polychromatic (even 
if possibly enhanced for k = ky) above the frequency of 
separation between the two ISB bare transitions, with 
angular separation of colors in free space. 

Of course, the present model is valid when the typical 
lifetimes of the cavity mode and of the electronic excita- 
tions are large enough to treat the different transitions 
independently. For what concerns the latter one, calcu- 
lations performed in a similar systeni^i show that the 
shortest lifetime (due to inelastic optical-phonon scatter- 
ing) is larger than 1.5 ps, corresponding to a broaden- 
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ing of about 0.5 meV. The cavity mode lifetime can be 
tuned by tailoring the optical cavity; in typical experi- 
ments the cavity mode broadening is of the same order 
of the electronic excitations one»i The different states in 
the intermediate and high-field regions are thus expected 
to be distinguishable in the experiments. Additionally, it 
is also worth recalling that in a typical reflectance spec- 
troscopy measurement, photons with fixed (i?, fc||) prop- 
agate in the substrate (of index Usuh) with fixed angle 9 
(with respect to the layers normal) given by: 



We show in Fig. [5] the energy versus plot for the 
polaritonic states around the resonance region of Fig. 0] 
(j^sub = 3.51, as in the experiments of Ref.|3j). The orange 
arrows connect the same states as in Fig.|4l i.e., those for 
which the central and generated polaritons have the same 
group velocities. As we can see, the angle difference is 
small (to ensure they all fall in the experimental light 
cone)^ but sizeable (slightly less than 10 degrees), so 
that all three states, even though broadened, couple to 
an external mode and can thus be in principle revealed 
in an experiment. 

The efficiency of the aforementioned non-linear pro- 
cess relies on inter-polariton interactions. In the context 
of exciton polaritons, polariton-polariton couplings have 
been previously studiedi^^— The corresponding scat- 
tering matrix elements depend however on the peculiar 



properties of the exciton components of the polaritons. 
The study of the scattering processes for intersubband 
polaritons in the presence of non-parabolic dispersions is 
however beyond the scope of this work. 



III. CONCLUSION 

In conclusion, we have shown that for intersubband po- 
laritons in narrow-gap semiconductors, with a significant 
non-parabolicity, the magnetic field plays a true role of an 
external control parameter that allows to tune the regime 
of light -matter interactions. It becomes then possible to 
tune the strength of the coupling of the light with the dif- 
ferent non-degenerate intersubband levels. We have re- 
ported numerical results for a InAs/AlSb system, and we 
have identified three different regimes for the polariton 
coupling function of the intensity of the magnetic 
field. Finally, we have presented a design for an opti- 
cal parametric oscillator in the FIR spectral range. The 
structure is based on the existence of a mirror dispersion 
scheme for the magneto-polaritons, which ideally allows 
fulfilling phase-matching requirements for the pump and 
parametric waves. 
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